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ABSTRACT 

We present the results of monitoring the radio emission from the Type lib 
supernova SN 200 Igd between 2002 February 8 and 2002 October 28. Most 
of the data were obtained using the Very Large Array at the five wavelengths 
of AA1.3 cm (22.4 GHz), 2.0 cm (14.9 GHz), 3.6 cm (8.44 GHz), 6.2 cm (4.86 
GHz), and 21 cm (1.4 GHz). Observations were also made with Giant Meterwave 
Radio Telescope at A21 cm (1.4 GHz). The object was discovered optically well 
after maximum light, making any determination of the early radio evolution 
difficult. However, subsequent observations indicate that the radio emission has 
evolved regularly in both time and frequency and is well described by the SN 
shock/circumstellar medium interaction model. 

Subject headings: galaxies: individual (NGC 5033) - radio continuum: stars - 
stars: mass-loss - supernovae: general - supernovae: individual (SN 200 Igd) 

1. Introduction 

Supernova (SN) 2001gd in NGC 5033 was discovered at magnitude, 16™5 on 2001 
November 24.820 UT by Nakano et al. (2001) at RA(J2000) = 13''13™23'!89, Dec(J2000) = 
-|-36°38'17'.'7, with independent discoveries by other observers following quickly. No source 
had been detected at that position (limiting mag. ~ 18.5) on a previous image taken on 2001 
May 8 and, nothing is visible on the Palomar Sky Survey red and blue plates (Nakano et al. 
2001). Ten days later, Matheson et al. (2001) confirmed the discovery as a type lib super- 
nova well past maximum fight, stating "the spectrum is almost identical to one of SN 1993 J 
obtained on day 93 after explosion." Using this suggestion to constrain the explosion date for 
the source, the optical discovery occurred 83 days after the explosion (Nakano et al. 2001) 
and; the first radio observation was 158 days after the explosion. Tully (1988) determined a 
distance of 18.7 Mpc (assuming Hq = 75 km s~^ Mpc~^) to NGC 5033. Accepting the Tully 
(1988) distance and using Hq = 65 km s~^ Mpc^^, we will adopt the distance to the SN of 
21.6 Mpc, which results in an absolute visible magnitude at discovery of < —15, although it 
is not known what the magnitude was at maximum optical brightness (Nakano et al. 2001). 
We use the Tully distance method to maintain uniformity with our distance determinations 
to other radio SNe including SN 1985L (Van Dyk et al. 1998). 

^National Research Council Postdoctoral Fellow. 

^On assignment from the Astrophysics Division, Space Science Department of ESA. 
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Stockdale et al. (2002) detected radio emission from SN 2001gd on 2002 February 8.54 
with the Very Large Array (VLA)^ at 22.485 GHz (A 1.3 cm), 8.435 GHz (A 3.6 cm), 4.885 
GHz (A 6.1 cm), and 1.465 GHz (A 21 cm). The radio position of RA(J2000) = 13'*13'"23!90, 
Dec(J2000) = +36°38'18'.'l, with an uncertainty of ±^.'2 in each coordinate (Stockdalc et al. 
2002), is in agreement with the optical position to within the errors. Since this first detection, 
we have made, and will continue to make, radio observations with the VLA. We note that 
NGC 5033 is also the host galaxy of two other historical supernovae, SNe 1950C and 1985L 
(Kowal, Sargent, & Zwicky 1970; Aksenov et al. 1985). The latter was also detected in the 
radio in 1986 (Van Dyk et al. 1998). Stockdale et al. (2002) did not detect radio emission 
from either of these historical supernovae, nor do we in any of our subsequent observations. 



The VLA radio observations reported here were made between 2002 February 8 and 2002 
October 28. VLA phase and flux density calibration and data reduction followed standard 
procedures (Weiler et al. 1986; Weiler, Panagia, & Sramek 1990), using 3C 286 as the 
primary flux density calibrator and J1310+323 as the phase calibrator and as a secondary 
flux density calibrator. 

In Table 1, we hst the flux density observations of the secondary cahbrator source 
J1310+323, which are plotted in Figure 1. The defined position for J1310+323 for phase 
reference is RA(J2000) = 13'*10"*28!66, Dec(J2000) = +32°20'43'.'8. 

In Table 2, we list the flux density measurements for SN 2001gd and Figure 2 plots the 
data. The data appear to be well described by a parameterized model described below and 
illustrated in Figure 2 . 

The flux density measurement error given for the VLA measurements in Table 2 is a 
combination of the rms map error, which measures the contribution of small unresolved 
fluctuations in the background emission and random map fluctuations due to receiver noise, 
and a basic fractional error, e, included to account for the normal inaccuracy of VLA flux 
density calibration and possible deviations of the primary calibrator from an absolute flux 
density scale. The final errors (cr/) as listed in Table 2 are taken as 



2. 



Radio Observations 




(1) 



^The VLA telescope of the National Radio Astronomy Observatory is operated by Associated Universities, 
Inc. under a cooperative agreement with the National Science Foundation. 
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where 5*0 is the measured flux density, ctq is the map rms for each observation, and e = 0.10 
at 21 cm, 0.05 at 6 cm, 0.05 at 3.6 cm, 0.075 at 2 cm, and 0.10 at 1.2 cm. 



3. Discussion 

Common properties of radio SNe (RSNe) include non-thermal synchrotron emission 
with high brightness temperature, turn-on delay at longer wavelengths, power-law decline 
after maximum with index /3, and spectral index a asymptotically decreasing to an optically 
thin value (Weiler et al. 1986). Weilcr et al. (1986, 1990) have shown that the "mini- 
shell" model of Chevalier (1982a; 1982b), adequately describes previously known RSNe. In 
this model, the rclativistic electrons and enhanced magnetic fields necessary for synchrotron 
emission are generated by the SN shock interacting with a relatively high-density envelope 
of circumstellar material (CSM) which has been ionized and heated by the initial UV/X-ray 
flash. This dense cocoon is presumed to have been established by a constant mass-loss (M) 
rate, in a constant velocity wind (unwind; i- ©■) Pcsm — ^^-^) from a red supergiant (RSG) 

'^^wind 

SN progenitor or companion. This ionized CSM is also the source of the initial absorption. 
The rapid rise in radio flux density results from the shock overtaking progressively more of 
the wind matter, leaving less of it along the line-of-sight to the observer to absorb the more 
slowly decreasing synchrotron emission from the shock region. 

3.1. Parameterized Model 

Following Weiler et al. (1986; 1990; 2001; 2002), Montes, Weiler, & Panagia (1997), 
and Chevaher (1998), we adopt the parameterized model. 



with 
where 



and 
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with Ki, K2, K2, and K4 corresponding, formally, to the flux density (i^i), homogeneous 
absorption {K2, K4), and clumpy (or filamentary) absorption at 5 GHz one day after 
the explosion date tn. The terms tcsm,, and tcsm , describe the attenuation of a 

■t ^ '--^'-'-'^^-'^homogcncous '--^'-'-'^^-'^ clumps 

local homogeneous and clumpy GSM that arc near enough to the supernova progenitor that 
they are affected by the expanding blastwave. The TcsMhomogoncous absorption term is due to an 
ionized medium which uniformly covers the emitting source. The TcsMdumps absorption term 
describes the attenuation produced by an inhomogeneous medium (see Natta & Panagia 
(1984) for a more detailed discussion). The Tdistant term describes the attenuation produced 
by a homogeneous medium that uniformly covers the source but is so far from the SN 
progenitor that it is not affected by the blastwave and is constant in time. 

All external and clumpy absorbing media are assumed to be purely thermal, singly 
ionized gas that absorbs via free-free (f-f) transitions with frequency dependence in 
the radio. The parameter 6 describes the time dependence of the optical depth for the 
local, uniform medium. For an undeceleratcd SN shock, S — —3 is appropriate^, and we 
have assumed this value for our parameterization since the data are insufficient to properly 
determine 6 at this time. The parameter 6' describes the time dependence of the optical 
depth for the non-uniform (clumpy) absorbing medium. For this source and at this epoch, 
we neglect the {K4) term because any distant, time-independent absorption effects will likely 
be evident only at lower frequencies and later times than are available in our current data 
set (see e.g. Montes et al. 1997; Schlegel et al. 1999 for SN 1978K). 

Ghcvalier (1998) proposed that synchrotron self-absorption (SSA) may be a significant 
source of absorption at high frequencies and early times for some RSNe. Given the relatively 
poor sampling of the radio emission turn-on at early times, SSA is not expected to be 
significant and our parameterization of SN 2001gd is not improved by the inclusion of SSA 
component. 

A search for the best parameter fit to Equations 2 through 6 was carried out by mini- 
mizing the reduced chi-squared (x^ed)- '^^^ best- fit parameter values and comparisons with 
the type lib SN 1993J and the other recent RSN in NGC 5033, SN 1985L, are fisted in Table 
3. The resulting model curves are plotted in Figure 2. The errors listed for the parameter 
values, determined by the "bootstrap" method, are at the la uncertainty level. 

The fitting procedure and examination of Figure 2 indicates that the gross properties 
of the radio emission from SN 2001gd are relatively well described by a purely thermal 
absorption model such as that given by Eqs. 2 — 6, with parameter values fisted in Table 3. 



"^Weiler et al. (1986) show that the absorption term TcsMh„mogonoous °^ ^ ^ ■ Since r oc for an undeceler- 
atcd shock, TcsMhomogeneous °^ ^'^^ from Equatlon 4, it follows that (5 = — 3 in this case. 
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The requirement for non-zero K2 and K3 parameter values implies the existence of a clumpy 
or filamentary absorption component embedded in a relatively homogeneous medium (on a 
more detailed level, this is consistent with Case 2 discussed by Weiler et al. 2001, 2002). The 
distant absorption term (K4) was set to zero for reasons discussed previously. The parameter 
values are typical of values determined for other type II SNe, although they are not tightly 
constrained by the relatively sparse data set. 

3.2. Comparison to Other RSNe 

It is useful to make distance independent comparisons between SN 200 Igd to SN 1985L, 
also in NGC 5033. Van Dyk et al. (1998) reported only two 6 cm detections of SN 1985L 
approximately 1 year after explosion, with numerous non-dctcctions. SN 1985L was identified 
as type IIL by Filippenko & Sargent (1985) and Kris (1985). Despite the relatively small 
amount of radio data available for the SN 1985L, there are some vahd comparisons that 
can be made. SN 2001gd (Lg cm peak = 2.91 x 10^''' erg s~^ Hz~^) is almost ten times more 
luminous at the 6 cm peak flux density than SN 1985L [Lq cm peak = 3.10 x 10^^ erg s~^ 
Hz~^; Van Dyk et al. 1998).^ The 6 cm peak luminosity of SN 200 Igd also lies within a 
range of other type II SNe, for example, SN 1979C (2.55 x 10^^ erg s'^ Hz-^) and SN 1980K 
(1.18 x 10^^ erg s~^ Hz~^). SNe 2001gd and 1985L also have very different times from 
explosion to 6 cm peak, although in both cases there is great uncertainty in this time scale. 
SN 2001gd has a much steeper spectral index [a = —1.4) than any other well-observed RSN. 
Weiler et al. (2002) demonstrate that most type II SNe have a spectral index {a > —1.0), 
similar to SN 1985L which has a moderate value oi a — —0.58. The power law dechne index 
for SN 2001gd (/3 0.96) is comparable to that of other type II RSNe. 

It is also useful to make some comparisons between SNe 200 Igd and 1993 J since both 
are type lib SNe and SN 1993J was much better sampled than SN 2001gd. Both have 
relatively equivalent 6 cm peak luminosities, as reported in Table 3. The time to 6 cm peak 
for SN 1993 J is in relatively close agreement with SN 2001gd (~ 170 days), although this 
value is somewhat atypical when compared with other type II SNe. Type Ib/Ic SNe typically 
achieve their 6 cm peak luminosity in <40 days, while most type II SNe take from as a few 
as ~100 days to well over 1 year (Weiler et al. 2002). Both SNe 2001gd and 1993 J have 
spectral indices < —1.0, which is more typical of type Ib/Ic SNe (Weiler et al. 2002). But, 



^We have adjusted the luminosity of SN 1985L to account for a Tully (1988) distance to NGC 5033 of 
21.6, assuming Hq = 65 km s'^ Mpc'S Van Dyk et al. (1998) assumed Hq = 70 km s"^ Mpc"^ for a Tully 
(1988) distance of 20 Mpc. 
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SN 2001gd appears to have a power-law decline index comparable to that of SN 1993J and 
other type II SNe, unlike values reported for type lb SNe 1983N and 1984L {f3 ~ —1.6; 
Sramek, Panagia, & Weiler 1984; Weiler et al. 1986; Panagia, Sramek, &Weiler 1986). 



3.3. Estimates of Mass-loss Rates 

Using the derived parameter values from the radio light curves, we can estimate the 
presupernova mass loss rate of the SN 2001gd progenitor star (Weiler, Panagia, & Montes, 
2001; Weiler et al. 1986; 2001; 2002), as 



M/(M0 yr-i) ^ ^ ,,,, / u, 



\ eff / g-i; 



1.5 



Wwind/(10 km s-1) ^ ^ VlO^ km s'^ 

UsdaysJ [tj [iO^k' ' 



where 



(^'efr'') - 0-67[(rcSMhomogeneou. + ^CSM.i^^p J ^ " '^CSMhomogeneou J '^CSMdumps 

This estimate depends on the parameterized model as well as on assumptions of tem- 
perature (TcsM = 2 X 10"^ K), pre-SN wind speed (t^wmd = 10 km s~^), shock velocity 
{vi = 10^ km s~^), and tj which accounts for the time between the appearance of optical 
lines from which the expansion velocity is determined, and the explosion date of the SN, 
typically assumed to be 45 days (Weiler et al. 1986). Equation 7 is derived from Case 2^ 
presented by Weiler et al. (2001; 2002). Since the SN was detected at such a late time, we 
can only approximately estimate the shock velocity and temperature at early times, making 
the mass loss rate estimates for SN 2001gd rather uncertain. However, the estimated mass 
loss rate of 3 x 10~^ M© yr~^ is close to the values found for both SNe 1993J and 1985L. 
As one might expect since they are both type lib SNe, there are many similarities between 
SNe 200 Igd and 1993 J, especially in that they require the presence of both uniform and 
clumpy (or filamentary) media to account for the radio light curve behavior while the mass 



^Case 1 is a uniform, external absorbing medium; Case 2 is a clumpy or filamentary medium with a 
statistically large number of clumps; and Case 3 is similar to Case 2, but with a statistically small number 
of clumps. 
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loss rate for SN 2001gd appears to be in the lower range of known mass loss rates for type 
II SNe (1 X 10"""' — 1 X 10~^ Mq yr"^), it clearly exceeds the values for typical Ib/Ic SNe 
(7 X 10-^—1 X 10-5 Mq yr-i; Weiler et al. 2002). 

4. Conclusions 

We present new observations of SN 2001gd at multiple radio wavelengths for the first year 
after explosion. Modehng implies that the evolution of the radio emission from SN 2001gd 

behaves in a systematic fashion and is consistent with the non-thermal emitting, thermal 
absorbing model previously used successfully for other RSNe. The CSM required by the 
model fits is best described as a "filamentary" or "clumpy" medium embedded in a uniform 
component and appears similar to that suggested for the type lib SN 1993J (Van Dyk et 
al. 1994). In general, SN 2001gd is very similar to SN 1993 J suggesting that type lib SNe 
may be relatively uniform in their properties. By comparison, the type Il-L SN 1985L in the 
same galaxy as SN 2001gd appears significantly different in its radio properties, illustrating 
the general diversity of type II RSNe. 

We are indebted to Dr. Barry Clark at the VLA for scheduling our numerous obser- 
vations and to observers who contributed their telescope time. KWW wishes to thank the 
Office of Naval Research (ONR) for the 6.1 funding supporting this research. CJS is a Na- 
tional Research Council Postdoctoral Fellow. We would like to thank A. Ray and P. Chandra 
for sharing their 21 cm GMRT observation. Additional information and data on radio su- 
pernovae can be found on http://rsd-unmv.nrl.navy.rnil/721Ji-/weiler/ and linked pages. We 
also would like to thank our anonymous referee for his/her useful comments. 
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Table 1. Calibrator Measurements for J1310+323 



Obs. 


Ref. Epoch 


Tel. 


S(21cm) 


S(6cm) 


S(3.6cm) 


S(2cm) 


S(1.2cm) 


Date 


(days) 


Config. 


(Jy) 


(Jy) 


(Jy) 


(Jy) 


(Jy) 


2001-NOV-24 


= 0.00 














2002-Feb-08 


75.72 


VLA-A 


1.703 


1.659 


1.765 




2.145 


2002-Mar-02 


97.60 


VLA-A 


1.703 


1.687 


1.745 


2.027 


2.164 


2002-Mar-14 


109.66 


VLA-A 


1.700 










2002-Mar-23 


118.42 


VLA-A 


1.681 




1.798 


2.077 


2.315 


2002-Apr-06 


132.38 


VLA-A 


1.690 


1.700*^ 


1.833*^ 


2.073 




2002-Jun-12 


199.30 


VLA-B 


1.668 


1.726 


2.005 






2002-Jul-16 


234.17 


VLA-B 


1.722 


1.773 


2.045 






2002-Oct-28 


337.15 


VLA-C 






2.130 


2.551 


2.728 



"^Interpolated from other calibrator measurements. 
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Table 2. Flux Density Measurements for SN 2001gd 



Obs. 


Rcf. Epoch 


Tel. 


S(21cm) ± (Tf 


S(6cin) ± o-j 


S(3.6cm) ± af 


S(2cm) ± cr/ 


S(1.2cm) ± <7f 


Date 


(days) 


Config. 


(mjy) 


(mJy) 


(mJy) 


(mJy) 


(mJy) 


2001-NOV-24'" 


= 0.00 














2002-Feb-08 


75.72 


VLA-A 


0.848 ± 0.129 


5.050 d= 0.257 


4.090 ± 0.208 




1.370 ± 0.176 


2002-Mar-02 


97.60 


VLA-A 


1.414 ± 0.177 


5.414 ± 0.284 


3.732 ± 0.195 


2.183 ± 0.242 


1.335 ± 0.182 


2002-Mar-14 


109.66 


VLA-A 


1.840 ± 0.352 










2002-Mar-23 


118.42 


VLA-A 


1.830 ± 0.339 




3.660 ± 0.208 


2.330 ± 0.322 


1.000 ± 0.187 


2002-Apr-06 


132.38 


VLA-A 


1.990 ± 0.344 


5.330 ± 0.277 


3.780 ± 0.212 


2.660 ± 0.304*^ 




2002-Jun-12 


199.30 


VLA-B 


3.230 ± 0.365 


4.842 ± 0.249 


3.100 ± 0.168 






2002-Jul-16 


234.17 


VLA-B 


2.733 ± 0.314 


4.260 ± 0.234 


2.390 ± 0.137 






2002-Sep-22 


301.40 


GMRT 


3.550 ± 0.220^^ 










2002-Oct-28 


337.15 


VLA-C 






2.071 ± 0.136 


0.893 ± 0.178 


0.809 ± 0.155 



^Nakano et al. (2001); Matheson et al. (2001); SN 2001gd is believed to have exploded roughly 83 days prior to its initial optical 

detection. See Section 1 for further explanation. 

''This value was not included in the modeling. 
'^Chandra, Ray, & Bhatnagar (2002) 



Table 3. Fit of SN 2001gd & Comparison with other Type II Supernovae. 



Parameter** 


SN 2001gd 


SN 1993J'' 


SN 19851^^ 




(1.49;^-^^) X 10^ 


4.14 X 10^ 


4.8 X 10^ 


a 


-1 qO + 0.13 


-0.99 


> -0.58 






-0.64 


< -0.71 


K2 


(3.25tt;i) X 


1.35 X 10^ 


> 7.00 X 10^ 


S 


= -3.00 


-1.99 


= -3.00 


K3 




2.54 X 10* 


0.00 


5' 


1 97+0.40 
J^-^' -0.15 


-2.02 


0.00 


Time to Lq cm peak (days) 
Le cm peak (erg Hz-^) 


17'Stlt 


168 


= 30 


(2.9ll°;J^) X 10^7 


1.45 X 10^^ 


~ 3.1 X 10^6 


il (M: yi;--' ) 


:-!.() X 10-" 


2.1 X 10-"' 


> G.:-') X 10-^ 



Assumes a constant mass-loss rate, constant velocity, wind-established CSM, 
i.e., p (X r~^. 

i-Van Dyk et al. (1994) 
'^Van Dyk et al. (1998) 
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Fig. 1. — Plot of the flux denstity of the VLA calibrator source J1310+323 at 1.3 cm (circles) , 
2 cm {crosses), 3.6 cm [squares), 6 cm (triangles), and 21 cm [diamonds). Note that the 
values for 21 cm are all reduced by 0.5 Jy to avoid confusion on the plot. 

Fig. 2. — The radio measurements and parameterized model for SN 2001gd at 1.3 cm [circles; 
solid line), 2 cm [crosses; dashed line), 3.6 cm [squares; dash-dot line), 6 cm [triangles; 
dotted line), and 21 cm [diamonds; dash-triple dotted line). The lines represent the best 
fit model as described in the text with the parameters listed in Table 2. SN 200 Igd is 
believed to have exploded 83 days prior to its initial optical discovery, see Section 1 for 
further explanation(Nakano et al. 2001; Matheson et al. 2001). 
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